Atomistic modelling of large-scale metal film growth fronts 
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We present simulations of metallization morphologies under ionized sputter deposition condi- 
tions, obtained by a new theoretical approach. By means of molecular dynamics simulations using 
a carefully designed interaction potential, we analyze the surface adsorption, reflection, and etching 
reactions taking place during Al physical vapor deposition, and calculate their relative probability. 
These probabilities are then employed in a feature-scale cellular-automaton simulator, which pro- 
duces calculated film morphologies in excellent agreement with scanning-electron-microscopy data 
on ionized sputter deposition. 

PACS: 81.10.Aj, 81.15.-z, 68.55.-a 
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In this Letter we present results of a hierachy of theo- 
retical models developed to describe the growth of metal 
thin films. Atomic-scale molecular dynamics (MD) and a 
feature-scale cellular-automaton simulator are combined 
to yield realistic simulations of film growth during physi- 
cal vapor metallization of contact vias typical of semicon- 
ductor device technology. The MD simulations, account- 
ing in full for the microscopic many-atom dynamics, are 
used to predict the reaction rates of the processes rele- 
vant in physical vapor deposition; the cellular automaton 
simulator incorporates the reaction rates thus obtained, 
enabling us to predict and understand the topography of 
/Ltm-scale film fronts, and their relationship to substrate 
geometry, incident beam energy, and angular beam dis- 
tribution. Our approach yields a consistent and compu- 
tationally efficient scheme to predict the topography of 
metal films on arbitrarily-shaped substrates. A compar- 
ison with scanning-electron-microscopy data on sputter- 
deposited Al-covered trenches demonstrates an excellent 
level of agreement between theory and experiment. 

Our theoretical approach proceeds in three steps, 
namely (a) a classical interatomic interaction potential 
for Al is developed, (b) reaction rates for Al atoms in- 
cident on Al surfaces are calculated therewith in a MD 
simulation, and (c) a cellular automaton is developed and 
employed to simulate /xm-scale film fronts, using the re- 
action rates extracted from MD. In the following, we an- 
alyze this procedure and present the relevant results for 
each step. 

Interatomic potential ^ -Previously developed classi- 
cal interatomic potentialaJTJ for Al-Al interactions have 
mostly focused on bulk and molecular properties. The 
Ercolessi- Adams potential! is a partial exception since, 
besides reproducing quantitatively the bulk and elastic 
properties of Al, it also yields interlayer relaxationsQ in 
good agreement with experiment E Using this potential, 
we obtain surface energies in excellent agreement with ab 
initio calculations.^ We thus chose this model as a start- 
ing point to develop a new classical many body potential 
for Al. Our own potential is carefully designed to repro- 
duce the properties of Al aggregates in a wide range of 



bonding configurations (from bulk Al to Al surfaces and 
Al molecules), with special regard to surface properties. 
Its functionalities significantly extend those of previous 
embedded atom models, while maintaining a high level of 
accuracy in reproducing Al bulk properties and surface 
formation energies.tj FixsL we introduced an exponential 
repulsive pair potentialB'tEl to account for the short-range 
interaction of Al atoms with kinetic energies exceeding 10 
eV; this is a key requirement, as the kinetic energies of de- 
posited atoms reach over 150 eV during ionized physical 
vapor deposition. Second, the embedding function has 
been readjusted to reproduce observed properties of low- 
density Al structures; after these changes, we obtained 
a significantly improved agreement of sever al , r eference 
quantitiestlin comparison to ab inito resultstHlEj and/or 
experiment O Third, we introduced a 5th-order polyno- 
mial cutoff function, smoothly cutting off the interaction 
range of our potential at an interatomic distance of 5.56 
A, slightly larger than third-nearest neighbor distance in 
bulk Al. 

As a further test of our potential, we calculated the 
barriers for homodiffusion on the low-index faces of fee 
Al. Specifically, we considered hopping on Al (110) along 
and orthogonal to the [110]-oriented rows, hopping on Al 
(100), and concerted exchange on Al (100). The gen- 
eral level of agreement with ab initio density-functional 
calculationsaiij is very good (as can be seen in Table 
||), giving us further confidence in the accuracy of our 
potential. 

TABLE I. Comparison of selected hopping and exchange 
diffusion barriers on low-index Al surfaces obtained with the 
present model and in ab initio calculations. Al (111) is also 
included for completeness. 



This work 



Ab initio 



Al (111) hopping 
Al (100) hopping 
Al (100) exchange 
Al (110) _L hopping 
Al (110) || hopping 



0.04 
0.60 
0.50 
1.13 
0.30 



0.04 a 
0.68 a ,0.65 b 
0.35 a 
1.06 a 
0.60 a 



^Reference |L 2l "Reference O 
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Reaction rates from molecular dynamics - In the sec- 
ond step of our approach, reaction probabilities were cal- 
culated in classical molecular dynamics simulations us- 
ing our Al interaction potential. In particular we deter- 
mined, as a function of the energy and off-normal angle 
of incident Al atoms, the probability of three processes: 
adsorption, reflection, and etching (in the latter, the in- 
coming atom's impact on the surface causes the kick-out 
of one or more substrate atoms). Supercells containing 
1320 atoms arranged in 10 atomic layers are employed; 
cell dimensions are chosen so as to avoid artifacts of the 
in-plane periodicity. The starting configuration is chosen 
to be a (111) surface, the one Al surface with the lowest 
formation energy. All atomic coordinates are allowed to 
evolve dynamically, except those of the two bottom lay- 
ers of the supercell. The surface temperature is set at 
450 K (i.e about 1/2 of the melting temperature, and 20 
% larger than the bulk Debye temperature). 

We start our simulations with the incident Al atom 
placed outside the interaction range of the surface. Its 
initial kinetic energy is set in the range of 25 to 125 eV, 
and its starting angle off the surface normal in the range 
0° to 60°, which corresponds to typical ionized physical 
vapor deposition conditions. The trajectories of the inci- 
dent atom, and of any other atom which may be etched 
away from the surface upon impact, are then monitored 
until either a certain time span has elapsed, or the out- 
coming atoms (in the case of reflection or etching) have 
traveled a distance of 10 A away from the surface. Ana- 
lyzing 200 trajectories per incident energy and angle, we 
collected a statistically significant sample of well-defined 
adsorption, reflection, and etching events. The relative 
probability of the corresponding process is calculated as 
the ratio of the number of events of each kind to the to- 
tal number. The typical statistical error in the reaction 
probabilities thus determined is below 5 %. 

The behavior of reflection, adsorption, and etch rates 
as a function of off-normal angle are summarized in Fig. 
|l| for two representative incident kinetic energies, namely 
25 eV in panel (a), and 100 eV in panel (b). At low en- 
ergies (panel (a)) the adsorption probability (thick solid 
line, solid squares) decreases from near unity for angles 
below 20° to approximately 1/2 for angles in the range 
of 60° . This decrease is compensated by a corresponding 
increase in the reflection probability (dashed line, solid 
triangles). At low energy, etch processes (thin solid line, 
solid circles) are negligible at all angles. Increasing the 
incident kinetic energy to 100 eV, we find considerable 
changes in the relative reaction probabilities. Even at 
small angles, the reflection probability is non- vanishing; 
the adsorption probability is correspondingly reduced, 
and is now below 0.7 at all angles. More interestingly, 
the etching probability is always non-zero, it reaches a 
maximum of 0.4 at 50°, then decreases as near-grazing 
angles are approached. For small deviations from normal 
incidence, the etch rate initially raises, since the proba- 
bility of a surface atom to gain momentum directed away 
from the surface increases when the incoming atom ar- 



rives at an oblique angle at the surface. At large angles 
of incidence the etch rate drops because of the competing 
specular reflection events. 
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FIG. 1. Reaction probabilities for Al atoms impinging on 
Al(lll) with a kinetic energy of 25 eV (panel (a)) and 100 
eV (panel (b)) as a function of the off-normal angle. The 
processes considered are adsorption (thick solid line, solid 
squares), reflection (dashed line, solid triangles), and etching 
(thin solid line, solid circles). 

By adjusting the bias voltage between sputter source 
and deposition target during ionized physical vapor de- 
position, both the energy regimes just discussed are ex- 
perimentally accessible. It is expected that they lead to 
rather different trench topographies, which we simulate 
with a cellular-automaton technique. 

Feature-scale cellular- automaton simulator - In the 
final step of our approach, we have developed a two- 
dimensional cellular automaton to model the growing 
film front on a /im scale. The automaton accounts for the 
effects of flux shadowing, adsorption, reflection, etching, 
and surface diffusion. The simulated structure is repre- 
sented in cross section by a two-dimensional grid. Each 
grid cell represents an Al atom, and is assumed to have 
a physical length of 2.5 A (the effective atomic diameter 
in Al bulk): thus, for example, a l-/i-wide structure will 
be described by 4000 grid cells across. 

The atoms are serially and independently emitted from 
the sputter source far above the surface according to a 
pre-determined angular and energetic distribution, and 
move on a straight trajectory, determined also by the 
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applied source-target bias, until they strike the growing 
film front. Interactions in the gas phase are neglected 
in view of the the low pressure (typically a few tenths 
of mTorr) and the resulting long mean free path typi- 
cal of sputter deposition. Spontaneous desorption is also 
negligible in all the conditions considered. 

The impact angle and energy of the atoms hitting 
the surface determine which process is activated upon 
impact.t^l This may be any one of the three (adsorption, 
reflection, etching) whose rates have been previously cal- 
culated via MD simulations. If the atoms are reflected, or 
an etch process takes place, the path of the corresponding 
atoms is further traced until they hit the film surface for 
a second time. The three basic processes can then take 
place again, and so forth. Finally, the atoms get either 
adsorbed, or escape back into the gas phase. 

Adsorption present some additional complications. 
The local diffusion and accomodation mechanism upon 
adsorption is a key ingredient in deposition models. A 
reasonable assumption is that incoming atoms are ac- 
comodated at a (local) minimum energy site within one 
diffusion length from the landing site, the diffusion length 
being determined by the surface temperature and mor- 
phology, and the deposition rate. In the cases of inter- 
est here, the problem of determining an effective diffu- 
sion length for the adsorbed atoms is quite formidable 
for several reasons. On the one hand, one expects low 
effective diffusivity due to the very high experimental 
growth rate (not well controlled, bill-in the order of 0.5 
/Ltm/minute, or roughly 40 ML/sect3'EZl) , and also because 
the growing surface rapidly looses its low-index charac- 
ter becoming essentially disordered. On the other hand, 
collision energies are large, and (although energy transfer 
mechanisms on rough surfaces are largely unknown) one 
may expect multi atom events and transient mobility ef- 
fects to increase the effective diffusivity. We aim at using 
parameters that maximize the coordination of the new 
particle, with the constraint that the impact-to-final site 
distance i&aninimized (so that the local film curvature is 
minimizedta) . 

In practice, we set the following criteria for (a) the 
maximum diffusion distance, and (b) the nature of a 
minimum-energy site: (a) The maximum diffusion dis- 
tance is taken to be d max = d^lx = 5 grid spacings =12 
A. All sites within this distance from the landing site are 
analyzed; the chosen final site is the one closest to the 
landing site, among those with the highest local coordi- 
nation (to be defined below, point (b)). Because of this, 
the effective diffusion length is actually rather smaller 
than the maximum value d max . Our specific choice of 
rfmax may be simply be regarded as a calibration factor 
designed to avoid dendritic structures or very flat sur- 
faces which are not observed in the experimental regime 
of interest here. However, we find a rather sharp change 
in the dependence of the surface roughness on d max , from 
rather weak for d max > d**^ to strong for d max < dj^.tj 
Thus our value d^ax i s effectively a crossover threshold 
between the two regimes mentioned. 



(6) - A minimum.-energy site is defined in terms of its 
local coordinationEJ as follows: for each candidate site 
within the preset maximum distance from the landing 
site (point (a)) we calculate the number of atoms con- 
tained in a circle center at the candidate site and with 
a radius of 7 grid spacings. This criterion amounts to 
selecting the site with maximum average coordination; 
while similar to that of highest first neighbor coordina- 
tion, contrary to the latter it avoids pathological choices 
such as high coordination sites on highly ramified struc- 
tures. 

Results of topography simulations - Figure ^ depicts 
trench topographies predicted by our model for different 
deposition conditions, compared to scanning-electron- 
microscope pictures taken in similar conditions.Ej The 
structure size is 1.2 /im across. In the cellular automa- 
ton simulation, the emitted atom energies are picked from 
a Thompson distribution centered at 3 eV, as suggested 
experimentallyOO The initial angles of the non-ionized 
atoms are chosen from a collimated cosine distribution 
with a maximum off-normal angle of 40° . To mimic the 
experimental conditions, we assume that 80% of the emit- 
ted Al atoms get one-fold ionized; for these atoms, the 
trajectory and impact energy change according to the 
applied source-target bias. 
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FIG. 2. Film morphologies on trench structures predicted 
for different ionized lpagnetron sputtering conditions com- 
pared to experiments. Il3 In panel (a) we set 80% ionization, 
and 10 V bias, in panel (b) 80% ionization, and 80 V bias. 
Below, a SEM micrograph with experimental results. 
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We now describe our results for different values of the 
the sputter source-to-target bias corresponding to the low 
and high energy regimes identified in our MD calcula- 
tions. The results in panel (a) of Fig. || were obtained 
setting the bias to 10 eV. In agreement with experiment 
(also displayed in Fig. 2), we predict a film growth front 
of rounded shape on top of the feature and, due to geo- 
metric shadowing, a reduced film thickness at the bottom 
of the trench. The pile-up at the center of the trench is 
not only of geometric origin, but also partly due to re- 
flections of atoms impinging on the trench sidewalls. In 
panel (b) of Fig. || we report results obtained with a bias 
of 80 eV. Our calculations predict, in accord with exper- 
iment, a roof-like structure on top of the feature. This 
structure is due to the preferential etching at angles of 
50° (see Fig. [l] (a)) which leads to a lower deposition 
rate on the roof-like structure. Simulations of structures 
scaled down in size by a similarity transformation)^ (i.e. 
with the same geometry, relative sizes, and aspect ratios) 
produced very close results, suggesting that the profiles 
obtained are largely self-similar at these length scales.EHl 

In conclusion, we have demonstrated the viability of 
accurate simulations of mesoscopic thin film morphol- 
ogy based on atomic-scale simulations. We performed 
detailed theoretical calculations of the probabilities for 
surface reactions taking place during ionized physical va- 
por deposition conditions, and combined these predic- 
tions with /xm-scale cellular-automaton simulations. Our 
molecular dynamics calculations revealed strongly energy 
dependent adsorption, reflection, and etch rates (the lat- 
ter exhibiting a distinct maximum for high incident ki- 
netic energies at ~ 50°). We were able to predict to- 
pographies of metal films deposited on trench structures 
under different ionized physical vapor deposition condi- 
tions in remarkable agreement with experiment. Our re- 
sults represent a major step ahead overr-sarlier thin film 
growth models based on rate equations^ which did not 
incorporate beam-energy-dependent surface reactions on 
an atomistic scale. 
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